Constrained simulations of the local Universe with modified gravity

We present a methodology for constructing modified gravity (MG) constrained simulations of the local Universe using positions and peculiar velocities from CosmicFlows data set. Our analysis focuses on following MG models: normal branch Dvali-Gabadadze-Porrati (nDGP) model Hu-Sawicki $f(R)$ model. develop model-independent with any given power spectra numerically calculated linear growth functions. Initial conditions (ICs) set are constructed standard cosmological $\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$ models. Differences between models' reconstructed Wiener filtered density resultant simulation presented showing importance generation ICs to study subtle effects in Universe. These first ever produced. The current work paves way improved approximate methods models scale-dependent functions, such as $f(R)$, high-resolution hydrodynamical zoom-in